Measurement-assisted Coherent Control 
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Two advantageous roles of the influence of measurement on a system subject to coherent control 
are exposed using a five-level model system. In particular, a continuous measurement of the popula- 
tion in a branch state in the Kobrak-Rice extended stimulated Raman adiabatic passage scheme is 
shown to provide a powerful means for controlling the population transfer branching ratio between 
two degenerate target states. It is demonstrated that a measurement with a large strength may 
be used to completely shut off the yield of one target state and that the same measurement with 
a weak strength can dramatically enhance the robustness of the controlled branching ratio against 
dephasing. 



I. INTRODUCTION 

Interest in coherent control of atomic and molecular 
processes 01 CJ [3- has grown rapidly in recent years. 
One of the prominent goals is the manipulation of quan- 
tum interference effects to enhance the yield of a par- 
ticular wanted final state while suppressing the yield of 
another. An example of this class of operations is en- 
hancement of the yield of one product of a reaction while 
suppressing the yield of a competing product. In this 
case, and in others as well, a quantum system under con- 
trol is typically not measured until the desired unitary 
evolution ends. The advantages and disadvantages of the 
influence of measurement on the system in a variety of 
coherent control scenarios remain unexplored. 

We note that the inhibition or acceleration of a quan- 
tum process under observation, the so-called quantum 
Zeno and anti-Zeno effects, demonstrate that the role of 
measurement in quantum dynamics can be dramatic. Al- 
though the quantum Zeno effect was first derived in terms 
of wavefunction collapse induced by a measurement, re- 
cent studies 0, B have shown that introducing the 
concept of wavefunction collapse is unnecessary for de- 
scribing the quantum Zeno and anti-Zeno effects. That is, 
a measurement apparatus, even in the case of obtaining 
a null result, always disturbs the system's Hamiltonian 
and can therefore strongly affect the quantum dynam- 
ics without involving our consciousness, a view that is 
adopted in this paper. 

In this paper we describe two interesting and advanta- 
geous roles of measurement in coherent control, using a 
well-studied model of selective photochemistry, namely, 
the Kobrak-Rice extended five-level stimulated Raman 
adiabatic passage scheme 0, 0, . That model is also 
the basis for the Chen-Shapiro-Brumer strong field con- 
trol approach to selective photochemistry Q ■ In particu- 
lar, it is shown that a continuous measurement of popu- 
lation in a so-called "branch state" with a large mea- 
surement strength provides a powerful means for ma- 
nipulating the nonadiabatic coupling and therefore the 
useful quantum interference between two particular adi- 
abatic states. The resultant control over the population 
transfer branching ratio between two degenerate product 



states is extraordinary insofar as (I) the yield of one of 
the two target states may be completely shut off while the 
yield of the other target state is still considerable, and (2) 
the intermediate state is unpopulated during the popula- 
tion transfer. Further, we show that the same measure- 
ment with a weak strength can dramatically enhance the 
stability of the controlled branching ratio under adverse 
situations, such as in the presence of strong dephasing 
modelled by stochastic energy level fluctuations. These 
results are of broad theoretical and experimental interest. 
Indeed, the measurement-assisted coherent control of the 
branching ratio between two unmeasured product states 
can be regarded as a third aspect of measurement effects 
that complements the quantum Zeno and anti-Zeno ef- 
fects. 




FIG. 1: A schematic diagram of the Kobrak-Rice five-level 
model [T(i|j for the control of population transfer branching ra- 
tio between two degenerate states |3) and |4), with the branch 
state |5) subject to a continuous population measurement. 
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II. EXPLOITING MEASUREMENT EFFECTS: 
STRONG MEASUREMENT LIMIT 

A schematic diagram of the Kobrak-Rice five-level 
model is shown in Fig. 2] The initial state |1) is cou- 
pled with the intermediate state |2) by a pump pulse, 
while state 2) is coupled with two degenerate prod- 
uct states |3) and |4) by a Stokes pulse. Due to the 
degeneracy of states |3) and |4), altering their popu- 
lation transfer branching ratio (denoted by B) cannot 
be achieved unless we introduce the "branch state" |5) 
that is also coupled with states |3) and |4) by a third 
"branching pulse" . All the laser fields are on resonance 
with the respective transitions |12| and are assumed to 
be Gaussian, with their pulse duration sufficiently short 
that, except for one case discussed below, the natural 
lifetimes of all the five states can be taken as infin- 
ity. To be specific, the associated Rabi frequencies (as- 
sumed to be real) [see Fig. ^ are chosen to be Sip = 
Slpexp[-(t-T) 2 /T 2 ], S1 S3 ,S4 = fW 4 ex P H 2 /T 2 ], and 
f2 B3 ,B4 = Slp 3 „B4exp[-0.5(i - 0.5T) 2 /T 2 ], where the 
pump and Stokes pulses have been counter-intuitively or- 
dered Q|. The characteristic magnitude of the various 
peak Rabi frequencies is represented by SI. The branch 
state is also subject to a continuous population measure- 
ment by coupling it to a continuum. Within the frame- 
work of the slowly varying continuum approximation |3j 
the disturbance of the system by the continuous mea- 
surement can be described by a finite and controllable 
lifetime of state |5) 0,0, with a decay rate constant (de- 
noted r) to be tuned via the strength of the system-probe 
interaction. Using the rotating wave approximation and 
the interaction representation, the system Hamiltonian 
is then given by H = H r + Hr, where H r denotes the 
on-resonance Hamiltonian in the absence of branch state 
measurement, 



(1) 



and (k\Hr\j) = —iTSksSjs is due to the continuous mea- 
surement. The eigenvalues and eigenvectors of H r , de- 
noted by Afc and |A*), k = 1 — 5, can be obtained ana- 
lytically. In particular, there exists one null eigenvalue 
(denoted Ai) with its eigenvector given by 

|ao = -i- [n 3B) o, -n P n Bil n P si B3 ,o] T , (2) 

iVi 

where Sl^p = (^S3^B4 — ^S4^S3)- Other eigenvalues of 
H r are given by 



eigenvectors are given by 
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where fc = 2 — 5, and Sl 2 f represents the sum of the 
squares of all the five Rabi frequencies. The associated 
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where Nk is the normalization factor. Note that |Ai) has 
a node on state |2), can fully correlate with the initial and 
final states, and can be far away from |Afc) {k = 2 — 5). 
In the absence of Hr, |Ai) is the key eigenstate for adi- 
abatic passage |4j without populating the intermediate 
state, and the associated population transfer branching 
ratio between the two degenerate target states |3) and 
|4) is given by B = Bx = Cl 2 B4 /Cl B3 . Interestingly, |Ai) 
[see Eq. (|5J)] also has zero overlap with the branch state, 
so it survives in the presence of Hr- However, the ex- 
act forms of the other eigenvalues and eigenvectors of H 
are too complicated to be useful. Below we use various 
perturbation theory techniques to reveal the interesting 
consequences for the system of introducing Hr ■ 

Here we consider the case of strong system-probe in- 
teraction, i.e., r >> SI. Then the system is dominated 
by the measurement and the laser-induced couplings can 
be regarded as a perturbation. In particular, H can be 
written as H = (H — V B ) + V B) where 

(k\V B \j) = S!p 3 (43<5j5 + SjaSks) + fis^jw^js + 5j4^*5)-(5) 

The eigenvalue-eigenvector structure of (H — V B ) can be 
easily obtained since state |5) therein is decoupled by con- 
struction. However, because (H — Vb) has two (degen- 
erate) null eigenvalues, a naive perturbation treatment 
is bound to fail. Our strategy is to first make a special 
superposition of the two null eigenstates of (H — Vb) to 
construct the eigenstate |Ai) that is known to exist in the 
presence of V B ■ Then second-order perturbation theory 
is applied to the subspace that is orthogonal to |Ai). We 
obtain three complex eigenvalues and one purely imagi- 
nary eigenvalue (denoted A2) that lies closest to Ai. As 
such, for the consideration of nonadiabatic effects below, 
only the coupling between |Ai) and the eigenvector asso- 
ciated with X' 2 (denoted | A2)) is important. Specifically, 



x' 2 = -i [si 2 P (si| 3 + si| 4 ) + n| B ] 2 /\r(Ntf 



S!p(S!s 4 S!p4 + SI53SIP3) 


^S3^S4^S4 — ^S3(^P + ^Si) 
^S3^S4^S3 — Slp4(Slp + Sl| 3 ) 

i [n P (n B3 + n B4 ) + si 2 B ) /r 



(6) 



(7) 



where N 2 is the normalization factor. Equation JfJJl shows 
that the peak value of i\' 2 is of the order of Sl 2 /r. Hence, 
in effect, the continuous measurement of the popula- 
tion of the branch state with a large strength modi- 
fies the spectrum of H r such that an eigenstate (|A' 2 )) 
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that lies close to |Ai) is created and the spacing between 
them (| A x — A' 2 |) becomes tunable. Equation © shows 
that | A2 ) also has a node on the intermediate state but 
does not overlap with the initial state at early times, 
and that the (final) B associated with |A 2 ) is given by 
B 2 = VL 2 m /Vl 2 Bi = X/B X . 

Let Ci(£) and C 2 (£) represent the quantum amplitudes 
on the adiabatic states |Ai) and | A2) , where £ = t/T. 
With all other eigenstates of H neglected the quantum 
dynamics reduces to that of a two-level system, 



dC 2 (0/dt = -i\' 2 TC 2 {£) 



(8) 
(9) 



where 02i(£) = (Ai |dA 2 /<i£) describes the nonadiabatic 
coupling between |Ai) and \\' 2 ). Using Eqs. © and © 
one finds 

OaxCO = -2n P nsB(^S3^B3 + n S i^B4)/(NiN^. (10) 

Consider then three different regimes of T. First, in 
the adiabatic limit, i.e., (FIT) 2 >> IT, the decay term 
— iX' 2 T in Eq. © dominates such that C 2 is always 
negligible and so is the population loss from |Ai). This 
gives B k, B\. Second, in the quantum Zeno limit, i.e., 
(tiT) 2 << IT, —iX' 2 T is negligible and |A X ) and |A' 2 ) 
become essentially degenerate, resulting in strong nona- 
diabatic effects that are beyond our control. Indeed, in 
this case the branch state is observed too vigorously: any 
transition to it is frozen and so it becomes useless for 
modulating control. The associated B then approaches 
^| 3 /f2| 4 , a result inherent to the four-level system ob- 
tained by discarding state |5). The third regime, i.e., 
((IT) 2 ~ FT, is of most significance. Denoting C\ and 
C 2 as the final values of Ci(£) and C 2 (£) obtained from 
Eqs. ©, ©, © an d (EH i an( i using Eqs. © and © 
one obtains 



B = 



BA 



2 I "S3 + (ClC| + C{C 2 )£lB?,£lBA t 



"S3 



c 2 \ 2 n 2 Bi - {dci + c^c 2 )n B3 n B4 



(11) 



Clearly, the value of B predicted by Eq. is not a 

simple mixture of the two branching ratios B\ and B 2 
associated with states |A) and |A 2 ). Rather, due to their 
constructive or destructive quantum interference B can 
be zero or infinity. That is, by varying F and therefore 
C\ and C 2 , different superpositions of the two adiabatic 
states |Ai) and |A' 2 ) can be realized and it becomes possi- 
ble to completely suppress population transfer to one of 
the two degenerate target states. Detailed conditions for 
B = or B = 00 will be discussed elsewhere jl3| . 



III. EXAMPLES 

Figure [21 shows two typical computational examples of 
(strong) measurement assisted control of B. As seen from 
Fig. [21 the theoretical results predicted by Eq. are 



in excellent agreement with those obtained by directly 
solving the Schrodinger equation. The effect of a contin- 
uous measurement of the population of the branch state 
is profound. In particular, in case (a) the yield of state 
1 3) (denoted P3) is essentially turned off for FT ~ 750 
whereas the yield of state |4) (denoted P4) is about 22%; 
in case (b) P4 is totally suppressed for FT ~ 900 whereas 
P3 is as large as 53%. From Fig. |2K it is also seen that 
even though state |4) is coupled more strongly with the 
measured branch state than state |3), P4 may, counter- 
intuitively, be less affected by the measurement than P 3 . 
We stress that the significant control shown in Fig. |3 
is achieved without populating the intermediate state. 
This is the case because the entire two-dimensional sub- 
space spanned by |Ai) and |A 2 ) has a node on state |2), 
and is effectively decoupled from other adiabatic states. 
Remarkably, our further calculations show that one may 
also achieve B « or B ss 00 even when the natural life- 
times of the product states are much shorter (e.g., T/5) 
than the pulse duration |l3j| . 
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FIG. 2: The final population in degenerate states |3) and 
1 4) (represented by P3 and P4, B — P3/P4) as a function of 
the measurement strength characterized by F. ClpT = 10, 
VlB'iT = 30, Q.B4.T = 50. In case (a) n S3 T = 30, £l S 4T = 70 
and in case (b) Qs3T = 60, fis4T = 40. Solid and dashed lines 
represent the theoretical results from Eq. 1111 . and discrete 
points represent the results obtained by directly solving the 
Schrodinger equation. 

As an aside we note that Hr can be thought of as 
an "imaginary detuning" from the branch state. Thus, 
the above formalism with slight modifications (e.g., with 
a real A 2 ) seems to suggest that the quantum interfer- 
ence between states |Ai) and |A 2 ) may also be manipu- 
lated if we detune the branching pulse by a very large 
amount. However, a large detuning often induces reso- 
nant or nearly resonant transitions to other states that 
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are not included in the five-level system. Moreover, due 
to the different nature of the associated nonadiabatic dy- 
namics, the complete suppression of one product channel 
based on detuning is found to be much less common (e.g., 
it does not occur in cases (a) or (b) in Fig. |2J), and if ob- 
tained has a far more limited tolerance to variation of the 
natural lifetimes of the product states than that achieved 
by a strong measurement. 



IV. WEAK MEASUREMENT LIMIT 

We now turn to the case of weak system-probe inter- 
action, i.e., r << il. Here Hr can be treated as a per- 
turbation to H r . The null eigenvalue Ai and the null 
eigenvector |Ai) still exist. Let the other eigenvalues and 
eigenvectors be A J! and \X' k ), k = 2 — 5. To the zeroth 
order of Hr one has |AJ!) = \Xk)- This yields, to the first 
order of Hr, 

»(A50 - 5R(A fe ); (12) 

3f(A' fc ') = -rxl(n B3 n S3 + n B4 n S4 ) 2 /Nl (13) 

The above perturbative treatment gives |A' fc ' — Ai| ~ 
l-^fc — Ai I , so the coherent population transfer can adi- 
abatically follow |Ai), B therefore equals B\, and the 
measurement at first glance seems to play no role. Nev- 
ertheless, the imaginary part of AJ! suggests that even 
if the undesired adiabatic states |AJ.') become populated 
due to some uncontrollable factors, such as the instabil- 
ity in laser phases or energy level fluctuations, they may 
still be absorbed away by the measurement which how- 
ever keeps the population in the null eigenstate intact. 
More significantly, we find a strong correlation between 
3(AJ!) and Db, where Db represents the difference be- 
tween B\ and the time-evolving branching ratio given by 
Using Eq. (@J one finds 

B (\ 2 k n S4 + n B3 n SB )2 n% 3 [ ' 

Comparing the above result with Eq. (|13() one sees that, 
in general, states with a larger \\ will give a worse 
B but will be absorbed by the weak measurement more 
quickly Interestingly, if £lsd,&B3 + &Si&Bi = 0, then 
Q(AJ,') = at all times for k = 2 — 5. But this presents 
no difficulty since Eq. 1)14(1 will give Db = as well. 
That is, in the cases where states |A'j!) are unaffected by 
the weak measurement of the population of the branch 
state, they necessarily give the right B. Clearly, then, 
a weak continuous measurement of the population of the 
branch state can serve as a convenient and effective "error 
correction" tool in realizing adiabatic passage. 



V. EXAMPLE: SELECTIVE POPULATION 
TRANSFER IN THE PRESENCE OF STRONG 
DEPHASING 

We now apply the finding of the previous section to 
selective population transfer in the presence of strong de- 
phasing [il,[iall3) where the dephasing (caused by, e.g., 
the coupling with other degrees of freedom, bimolecular 
collisions, or solvent effects on a solute molecule) is simu- 
lated by stochastic Gaussian fluctuations 8u>k (A; = 1 — 5) 
associated with each of the five levels. The degener- 
acy between states |3) and |4) is assumed to be main- 
tained, i.e., 8lo 3 = 5lo 4 . The mean properties of 5u>k 
are chosen to be (Suik) — and (5u)k{t)5uJk> (£')) = 
5kk' A 2 exp[— \t — t'\/r], k, k' = 1 — 3, 5. A computa- 
tional example is shown in Fig. |3J where AT = 15 (a 
value that represents extremely strong dephasing) and 
t = 0.02T (a value that is comparable to 1/Cl to fur- 
ther enhance the dephasing-induced nonadiabatic effects 
[T^.lr^l '). The goal in this example is to achieve B = 25.0 
and it is assumed that there indeed exists a branch state 
that meets this goal, i.e., B\ — 25.0. As is seen in Fig. 
Et, the value of B without the measurement is terribly 
degraded by dephasing and is as small as 2.6. However, 
upon introducing a weak measurement of the population 
of the branch state, we obtain a dramatic increase to 
B = 24.0 for TT = 3 (Fig. EJd), B = 24.9 for TT = 6 
(Fig. Ot) and for TT — 9 (Fig. [Ji). As seen in Fig. 
13 this almost perfect recovery of the desired value of B 
is attained with a considerable total yield of the prod- 
ucts and with a wide range of the (weak) measurement 
strength. We have obtained similar results in numerous 
other cases, such as those with different t, B\, ^53/^54, 
or different timing and duration of the branching pulse. 



VI. CONCLUDING REMARKS 

It is important to note that although the measurement 
of the population of the branch state absorbs away some 
population by, for example, breaking up some molecules 
from an ensemble of molecules that are subject to the 
control fields, our purpose in introducing it is to gener- 
ate desired quantum interference effects. In the strong 
measurement case considered in Sees. II and III, the 
measurement creates a second adiabatic state that can 
be significantly populated during the population trans- 
fer and thereby induces quantum interference effects that 
are absent in the unmeasured five-level STIRAP system. 
In the weak measurement case considered in Sees. IV 
and V, the measurement of the population of the branch 
state enhances the robustness of the quantum interfer- 
ence effects inherent in the unmeasured five-level STI- 
RAP system. That is, during the population transfer 
those molecules that are broken up by the weak mea- 
surement have the potential to generate a wrong branch- 
ing ratio, and those that survive the weak measurement 
carry the desired quantum interference induced by the 
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FIG. 3: Measurement-assisted recovery of a desired branch- 
ing ratio (jBi = 25.0) in the presence of strong dephasing. 
Shown here is the time dependence of the population in states 
|3) and |4) (represented by P3 and P4) without (case a) mea- 
surement, or with measurement for (b) FT — 3, (c) FT — 6 
, and (d) FT = 9. & P T = 20, f2 S3 T = 50, ft S 4T = 40, 
&B3T = 15, Qb4,T = 75. The results are obtained by av- 
eraging over 1000 realizations of the stochastic energy level 
fluctuations. 



control, with the price that the total yield of the product 
states is less than 100%. Nevertheless, one should never 
let the best be the enemy of the good. This price is ac- 
ceptable, particularly when the main objective of coher- 
ent control is to alter at will the branching ratio between 
two product states. The results of our study can also 
be regarded as an extension of our previous work |17| . 
where we have shown that the decay of a target state or 
measuring the population of the target state can greatly 
improve the performance of a three-level STIRAP system 
in the presence of strong dephasing 0, 0, 0] . 

In conclusion, we find that the population transfer 
branching ratio between degenerate product states may 
be significantly altered by observing the system. Our 
results suggest that coherent control can be even more 
powerful than previously thought if we combine quantum 
interference effects with measurement effects. The effects 
we have described depend on the character of the influ- 
ences induced when the population of the branch state is 
measured, and are different from the kind of final state 
branching ratio that develops as the result of incoherent 
kinetic competition between parallel pathways. Poten- 
tial applications of this work include control of molecu- 
lar chirality, control of photodissociation, and control of 
isomerization reactions in solution. 



strong control fields. 

Hence, while a measurement on a quantum system sub- 
ject to coherent control necessarily introduces decoher- 
ence to the system, it can still be beneficial for coherent 
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